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ABSTRACT 

We have used the radial velocity variations of two sdB stars previously reported to 
be binaries to estabHsh their orbital periods. They are PG 0940+068, (P=8.33d) and 
PG 1247+554 (P=0.599d). The minimum masses of the unseen companions, assuming 
a mass of 0.5Mq for the sdB stars, are 0.090 ± O.OGSMq for PG 1247+554 and 0.63 ± 
O.O2M0 for PG 0940+068. The nature of the companions is not constrained further 
by our data. 
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1 INTRODUCTION 

Sub-dwarf B (sdB) stars are thought to be core helium burn- 
ing stars with core masses of about O.SM© and extremely 
thin {<, 0.02Mq) hydrogen envelopes (Heber 1986; Saffer 
et al. 1994). It is thought that the eventual fate of an sdB 
star is to cool to form a white dwarf with a mass of about 
O.5M0, which is low compared to the typical mass for white 
dwarfs (Bergeron, Saffer & Liebert 1992). The formation 
of low mass white dwarfs, and by implication sdB stars, is 
thought to involve a common envelope phase, in which a 
companion to a red giant star in engulfed by the expanding 
outer layers. The resulting friction results in the companion 
spiraling in towards the core of the red giant, ejecting the 
envelope at the expense of orbital binding energy (Iben & 
Livio 1993). This provides a natural mechanism for the loss 
of the surface hydrogen layers in sdB stars. 

The binary fraction of sdB stars is expected to be high 
given the scenario outlined above. AUard et al. (1994) found 
that 31 of their sample of 100 sdB stars show flat spectral 
energy distributions which indicate the presence of compan- 
ions with spectral types in the range late-G and early-M. 
They infer a binary fraction for main-sequence companions 
of 54 -66 percent, although the companions in their sur- 
vey are over-luminous compared to normal main-sequence 
stars. A similar conclusion was reached by Jeffery & Pol- 
lacco (1998) based on the detection of spectral features due 
to cool companions. Companions to sdB stars can also be 
detected in eclipsing systems such as HW Vir (Wood & Saf- 
fer 1999), PG 1336-018 (Kilkenny et al. 1998) and the sdB- 
white dwarf binary KPD 0422+5421 (Koen et al. 1998). Saf- 
fer, Livio & Yungelson (1998) found that at least 7 of their 
sample of 46 sdB stars show radial velocity variations. Sev- 



eral of these binaries have subsequently had their orbital 
periods determined (Moran et al. 1999), although further 
observations are required to determine the nature of the 
companions in these binaries. 

In this paper we present two new orbital period deter- 
minations from radial velocity measurements for binary sdB 
stars from the survey of Saffer, Livio & Yungelson (1998). 



2 OBSERVATIONS AND REDUCTIONS 

Most of the data for this study comes from observations 
obtained with the 2.5m Isaac Newton Telescope (INT) on 
the Island of La Palma. Spectra were obtained with the in- 
termediate dispersion spectrograph (IDS) using the 500mm 
camera, a 1200 line/mm grating and a TEK CCD (charge 
coupled device) as a detector. 

Additional spectra of PG 0940+068 were obtained us- 
ing three other instruments. Ten spectra were obtained with 
the grating spectrograph on 1.9m telescope at the South 
African Astronomical Observatory (SAAO). Grating No. 5 
(1200 lines/mm) and a 1.2arcsec slit were used with the SITl 
CCD detector. Four spectra were obtained with RGO spec- 
trograph and 82cm camera on the Anglo-Australian tele- 
scope (A AT). Full details of these spectra can be found in 
Maxted & Marsh (1999). 

Two additional spectra of PG 1247+554 were obtained 
with the dual beam spectrograph ISIS on the 4.2m William 
Herschel Telescope (WHT) on the Island of La Palma. We 
used a 12001/mm grating on the red arm to obtain spectra 
of Ha and Hel 6678. The detector was a TEK CCD and the 
slit width was 1.2arcsec. 

For all these instruments and detectors, sensitivity vari- 
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Table 1. Summary of the spectrograph/telescope combinations 
used to obtain spectra for this study. The slit width used in each 
case is approximately larcsec. 



Date 



Telescope 



Spectre- Resolution Sampling 







graph 


(A) 


(A) 


Jun97 


INT 


IDS 


1.0 


0.4 


Feb 98 


INT 


IDS 


1.0 


0.4 


Mar 98 


SAAO 1.9m 


- 


0.8 


0.4 


May 98 


AAT 


UCLES 


0.10 


0.04 


Jun 98 


AAT 


RGO 


0.7 


0.3 


Feb 99 


INT 


IDS 


1.0 


0.4 


Jul 99 


WHT 


ISIS 


0.73 


0.4 



ations were removed using observations of a tungsten cali- 
bration lamp. Bias images show no signs of any structure 
in the region of the spectrum so a constant bias level de- 
termined from a clipped-mean value in the over-scan region 
was subtracted from all the images. The exposure times used 
varied from 300s to 1200s. 

Extraction of the spectra from the images was per- 
formed automatically using optimal extraction to maximize 
the signal-to- noise of the resulting spectra (Home 1986). Ev- 
ery spectrum was bracketed by observations of the internal 
arc lamp at the position of the star. The arcs associated 
with each stellar spectrum were extracted using the profile 
determined for the stellar image to avoid possible systematic 
errors due to tilted spectra. The wavelength scale was deter- 
mined from a polynomial fit to measured arc-line positions. 
Uncertainties on every data point calculated from photon 
statistics are rigorously propagated through every stage of 
the data reduction. 

Finally, a single spectrum of PG 0940-1-068 was obtained 
with the cross-dispersed echelle spectrograph UCLES on the 
AAT. Details of the reduction of these spectra are similar to 
those described in Gatti et al. (1997). 

Details of all these instruments, the dates of all the ob- 
serving runs and the resolution and sampling obtained are 
given in Table hi 



3 PG 0940+068 

This star was first noted as a potential binary by Bragaglia 
et al. (1990) and later confirmed as such by Saffer et al. 
(1998). Saffer et al. also noted the presence of the He I 6678 
absorption line in their spectra which identifies this star as 
an sdB star rather than a white dwarf. 

To measure the radial velocities we used least squares 
fitting of a model line profile. This model line profile is 
the summation of three Gaussian profiles with different 
widths and depths but with a common central position 
which varies between spectra. The resolution of each instru- 
ment is accounted for in the fitting process. Only data within 
2000 km/s of the Ha line is included it the fitting process. 
We first normalise the spectra using a linear fit to the con- 
tinuum either side of the Ha line. We used a least-squares 
fit to a single spectrum to establish the shape of an initial 
model line profile. A least squares fit of this profile to each 
spectrum in which the position of the line is the only free 
parameter yields an initial set of radial velocities. The peri- 
odogram of these velocities shows the best period to be near 
8.33d. Other peaks in the periodogram are much less signif- 



icant. We used a sine fit to the radial velocities to fix the 
position of the Ha line in each spectrum in a simultaneous 
least squares fit to all the spectra to establish an optimum 
model line profile. We then re-measured the radial velocities 
of each spectrum using this optimum model line profile. The 
quality of the model line fits to the spectra estimated by-eye 
and from the chi-square statistic were generally good. These 
are the radial velocities given in Table 0. 

The parameters of a circular orbit fit to the radial ve- 
locities are given in Table M. The value of chi-squared is 
rather high. This is, perhaps, not surprising given the diverse 
sources of spectra on which this orbit is based, but might 
also suggest that there are additional sources of uncertainty 
in our measured radial velocities, e.g., image motion within 
the slit, instrument fiexure and telluric absorption lines. We 
have adjusted the uncertainties in the measured parame- 
ters accordingly, i.e., multiplied them by •v/l41. 1/(41 — 4). 
The measured radial velocities and the circular orbit fit are 
shown in Fig. hi. 

We also established the parameters of the orbit by fit- 
ting all the spectra simultaneously using a least squares fit 
in which the position of the absorption line is determined by 
a circular orbit. The parameters of the orbit and the model 
line profile are all free parameters in the fitting process. This 
gives essentially the same result as fitting the individual ra- 
dial velocities. We also subtracted the model profiles from 
the observed spectra to look for spectral features due to a 
companion star, but without success. 



4 PG 1247+554 (GD 319) 

This star was first noted as a potential binary by Saffer et al. 
(1998), who also noted the strong He I 6678 line which sug- 
gests that this star is a blue horizontal branch (BHB) star or 
low-gravity sdB star. There is a K star companion 2.7arcsec 
away but this was shown to be unrelated to PG 1247+554 
by McAlister et al (1996), i.e., it is an optical double star. 

We acquired 11 spectra of this star and its nearby com- 
panion star over 4 nights of March 1999 covering the Ha line 
and the He I 6676 line, 4 spectra over two nights of June 1997 
covering the Ha line only and two spectra covering the Ha 
line and the He I 6676 line on a night of July 1999. 

We same procedure described above for PG 0940+068 
to establishing a model line profile and measure the ra- 
dial velocities from the Ha line. These are given in column 
2 of Table M. We also measured radial velocities from the 
HeI6678 line in our 1999 spectra using a similar procedure 
but using only two Gaussian profiles to model the line. The 
velocities are given in column 3 of Table 0. 

The parameters of a circular orbit fit to the radial veloc- 
ities measured from both spectral lines are given in Table H 
and the fit is shown in Fig. H. The value of chi-squared is, 
again, rather high so we have adjusted the uncertainties in 
the measured param eters accordingly, i.e., multiplied them 
by 1^/45.3/(30 — 4). A simultaneous fit to all the spectra to 
determine the orbit gives essentially the same results and 
there is no trace of a companion star in the residual spectra. 
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Table 2. Radial velocity measurements for PG 0940+068 



HJD 


Velocity 


HJD 


Velocity 


HJD 


Velocity 


HJD 


Velocity 


-2450000 


(kms-i) 


-2450000 


(kms-i) 


-2450000 


(kms-i) 


-2450000 


(kms-i) 


852.4349 


47.6±2.5 


855.4455 


-58. 2± 2.9 


882.3250 


-75.5± 7.7 


969.8409 


34.6± 4.1 


852.4463 


60.0±2.6 


855.4511 


-61. 8± 3.1 


882.3323 


-56.7± 8.1 


1238.4055 


-64.2± 7.1 


852.4589 


43.8±2.2 


855.4561 


-62. 9± 3.1 


882.3397 


-71.4± 8.4 


1240.3832 


-69.3± 2.0 


852.6349 


39.7±2.6 


855.4611 


-58. 3± 3.1 


882.3484 


-49.2± 9.2 


1240.4840 


-61. 9± 1.8 


852.6434 


34.9±3.4 


855.5102 


-61. 6± 3.2 


882.3567 


-59.2± 8.6 


1240.6804 


-56.5± 3.4 


853.5760 


11.2±2.2 


855.5151 


-62. 7± 3.1 


882.3640 


-57.8±10.8 


1241.3896 


-27.4± 2.1 


853.5847 


3.5±3.2 


881.3208 


-131.8±18.1 


941.8613 


-13.4± 1.6 


1241.5396 


-15.4± 1.9 


853.6816 


6.4±2.5 


882.2942 


-92. 5± 8.3 


968.8381 


41. 9± 3.4 


1241.6821 


-13.8± 3.2 


853.6912 


3.1±2.4 


882.3016 


-71. 5± 9.0 


968.8419 


43. 1± 3.2 


1242.3831 


17.0± 2.1 


854.4479 


-26.4±3.8 


882.3089 


-72. 7± 8.1 


969.8371 


36. 4± 3.4 


1242.6059 


20.5± 1.9 


854.4519 


-15.7±3.7 















Figure 1. Measured radial velocities for PG 0940+068 and the circular orbit fit. 
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Table 3. Circular orbit fits. The systemic velocity is denoted by 
7 and the semi-amplitude by K. 



Parameter 




PG 0940+068 


PG 1247+554 


HJD{To)-2450000 


1000.066 ± 0.027 


1240.855 ± 0.002 


Period(d) 




8.330 ± 0.003 


0.602740 
± 0.000006 


7 (kms-i) 




-16.7+ 1.4 


13.8 ± 0.6 


K (kms-i) 




61.2 + 1.4 


32.2+ 1.0 


Mass function(M0) 


0.20 + 0.01 


0.0021 + 0.0002 


X^ 




141.1 


45.3 


n 




41 


30 



5 DISCUSSION 

If we assume a canonical mass of O.5M0 for the visible 
components of these binaries (Saffer et al. 1994), we can 
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calculate a minimum mass of 0.63 ± O.O2M0 for the un- 
seen companions to PG 0940+068 and 0.036 ± O.OOSMq for 
PG 1247+554. If the unseen companion to PG 0940+068 
were a main-sequence star it would have an absolute magni- 
tude, Mv = 8.3±0.9 (Henry & McCarthy 1993), where the 
uncertainty reflects the error in the mass and the scatter 
in the mass - luminosity relation. On average, the absolute 
magnitude of an sdB star is Mv = 4.5 ± 0.8, so a main- 
sequence companion is likely to be contribute at least Iper- 
cent of the light at V. While this would not have been seen 
in our spectra, it should allow the nature of the companion 
to be established through high signal-to-noise spectroscopy, 
perhaps around the Ca II IR triplet as in Jeflery & PoUacco 
(1998). 
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Figure 2. Measured radial velocities for PG 1247+554 and the circular orbit fit. 
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Table 4. Radial velocity measurements for PG 1247+554 



HJD 


Ha 


HeI6678 


-2450000 


(kms-i 




622.4567 


17.7 ± 2.1 


— 


622.4652 


15.7 ± 2.4 


— 


623.4270 


-11.3 ± 2.8 


— 


623.4370 


-18.6 ± 2.8 


— 


1240.5436 


17.0 ± 2.3 


15.2 ± 2.0 


1240.5541 


16.0 ± 2.2 


13.0 ± 2.0 


1240.6456 


-13.2 ± 2.2 


-17.7 ± 1.9 


1240.7433 


-14.3 ± 2.3 


-12.9 ± 2.1 


1241.4975 


23.9 ± 2.4 


29.0 ± 2.1 


1241.7215 


25.7 ± 2.0 


27.4 ± 1.8 


1241.7788 


5.5 ± 2.3 


2.7 ± 2.0 


1242.5637 


-12.7 ± 2.1 


-14.1 ± 1.9 


1242.6808 


19.2 ± 2.2 


19.5 ± 2.0 


1242.7546 


39.7 ± 2.3 


43.0 ± 2.0 


1243.6586 


-9.4 ± 2.1 


-11.5 ± 1.9 


1368.3812 


3.1 ± 1.9 


-3.2 ± 1.6 


1368.3877 


0.7 ± 1.9 


-2.8 ± 1.6 



6 CONCLUSION 

We have measured the orbital periods for two sdB binary 
stars and calculated the minimum mass of the companion 
star. The orbital period of PG 0940+068 is 8.33d and the 
minimum mass of the companion is 0.63 ± O.O2M0. The or- 
bital period of PG 1247+554 is 0.598655d and the minimum 
mass of the companion is O.O87±O.OO3M0. If the companion 
to PG 0940+068 is a main sequence star, it will contribute 



Ipercent of the hght at V, which may make it directly de- 
tectable. 
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